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ABSTRACT 

Most Fermi GRB spectra appear as either a broken power law ex- 
tending to GeV energies or as a broken power with a separate GeV 
power law component. Here we show that such spectra can be un- 
derstood in terms of magnetically dominated relativistic jets where a 
dissipative photosphere produces the prompt MeV emission, which is 
extended into the GeV range by inverse Compton scattering in the ex- 
ternal shock, with possible contributions from a reverse shock as well. 
The bulk Lorentz factors required in these models are in the range 
of 300-600, and the MeV-GeV time delays arise naturally. In some 
cases an optical flash and a sub-dominant thermal component are also 
present. 



Introduction 



The GRB spectra observed with the LAT detector on Fermi reveal a diver- 
sity of appearances. Some have spectra consistent with a single canonical Band 



function (Band et al. 1993) extending to the highest detected multi-GeV energies, 



while in others at the higher end of the Band spectrum a second, flatter spectral 
component emerges, which sometimes appears to have a cutoff ( Pelassa|2011 Pe'er 



2011). The origin of such a dichotomy, as well as the extent of the applicability of 



such spectra for deriving lower limits on the bulk Lorentz factor (Kocevski & The 



Fermi Collaboration 2012) are the subject of ongoing research and discussions. 



Here we investigate possible scenarios for these two spectral types in the con- 
text of magnetically dominated, baryon loaded outflows, where the prompt MeV 
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emission arises in the jet photosphere. Magnetically dominated models have re- 



cently received increased attention 


Drenkhahn & Spruit||2002; 


Lyutikov & Bland- 


ford|2003 Giannios & Spruit 2007; 


Giannios 2011; Komissarov et al. 2009 |McK- 


inney & Uzdensky 2011| Metzger et al. 2011a; Meszaros & Rees 2011; Bosnjak 


& Kumar 


2012 


Zhang & Yan 


2011 


), while photospheric emission, either in bary- 



onic or magnetic models where dissipation occurs near the photosphere, have been 



shown to be able to produce Band-like non-thermal spectra ( 


Pe'er et al. 


2006 


Lazzati et al. 


2011; 


Beloborodov 


2010 


Ryde et al. 


2010 


Pe'er et al. 


2012 


). In 



the scenarios that we discuss here we combine these two elements and add a third 
one, assuming that the prompt spectrum is a combination of the emission from 
two distinct regions, namely, a magnetically dominated dissipative photosphere 
and the external shock. We argue that the photosphere is initially responsible for 
the low energy ( GBM ) part of the Band-like prompt spectrum, the final observed 
spectrum being modified by inverse Compton scattering in the external shock, 
which results in a delayed high energy {LAT ) component. The latter includes a 
forward and possibly also a reverse shock contribution, dominated by self-inverse 
Compton (SSC) as well as external inverse Compton (EIC) of the photospheric 
photons. We consider here mainly the radial structure of the outflow, and restrict 
ourselves to leptonic mechanisms. 

Our primary goal here is to explore the generic features of such models, and to 
test the extent to which they qualitatively lend themselves to the interpretation of 
LAT spectra. That is, we do not perform detailed fits to individual bursts, which at 
this point would be premature given the preliminary nature and the uncertainties 
in the models. The purpose is to explore the connection between elements of the 
physical model and the presence or absence of broad spectral features, and how 
these depend on typical burst parameters. We concentrate on the prompt emission, 
from the time of the trigger up to times comparable to the deceleration time, when 
the external forward shock (and the reverse, if present) have fully developed. The 
latter should be representative of the physics in the late afterglow, which we do not 
address specifically here. We then qualitatively compare our results to several of 
the archetypal LAT bursts, such as GRB 090902B ( Abdo fc the Fermi collaboration! 
2009), which has a Band component and an extra high energy power law (PL) 



second component; GRB 090926, which also has an extra spectral component and 
shows a high energy cutoff ( Ackermann fc the Fermi collaboration|20lT ); and GRB 



080916C, which shows a single Band component extending to multi-GeV (Abdo 
et al. 1120091). 
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We show that in the context of these models, this variety of appearances can 
be attributed to differences in the strength of the various radiation components 
under plausible variations of the input physical parameters. Previous investiga- 



tions (e.g., Toma et al. 2011b Fraija et al. 2012) have considered a connection 



between the extra high energy component and inverse Compton and/or a reverse 
shock in a baryon dominated outflow. Here the magnetically dominated dynamics, 
photospheric input spectrum and consideration of separate SSC and EIC compo- 
nents from both forward and reverse shock result in a different and broader range 
of possible outcomes. Depending on the parameters, the SSC or EIC from the for- 
ward shock and the reverse shock (if the latter is present) can either result in an 
extra high energy component, or sometimes in an apparent single Band component 
extending to the highest energies. We show that, even if there is a pair-production 
cutoff in the nonthermal prompt emission at a few or a few tens of MeV, the ex- 
ternal inverse Compton radiation from the reverse shock can constitute a natural 
extension of the absorbed emission to form a continuous Band spectrum out to 
higher energies. We argue that if a reverse shock does not develop, the domi- 
nant high-energy components will be the forward shock SSC or EIC, resulting in 
a distinct extra component. For some parameters, however, the latter can also 
be produced by a reverse shock. A further result from these calculations is that, 
since the high energy component arises mainly in the external shock, where the 
compactness parameter is low, these high energy photons do not constrain via 77 
cutoff considerations the bulk Lorentz factor to as high values as considered in 
many of the recent analysis papers. 

In £j2] we discuss the basic magnetic jet model used. In ^3] we give the details 
of the dissipative photospheric (prompt) spectrum. In £|4] we discuss the various 
external shock synchrotron and inverse Compton radiation components. In §[5] 
we present numerical results for various cases. In §[6] we discuss our results and 
conclusions. 



2. Magnetic Model 

We assume a total luminosity L t = 10 53 Lt i53 erg/s is released at a radius 
7o = 10 7 ro,7 cm, lasting for 10 s in the central engine frame . This is initially (at 
To) mainly magnetic, which at larger radii gets gradually converted into bulk kinetic 
energy and non-thermal radiation. We assume that at the dissipative photosphere 
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a fraction ( r of L t is released as prompt radiation, whose spectrum is discussed 
in the next section. A nominal value for this here is ( r ~ 0.5, which we assume 
to consist mainly of a non-thermal component C, n t ~ Cr giving the Band power 
law (PL) spectrum (a smaller fraction of ( r is in a thermal component (th < Cnt, 
discussed in £3.2, so that ( r = ( nt + (th)- The remaining energy fraction beyond 



the photosphere (k, after the Lorentz factor reaches its saturation (coasting) value, 
is in kinetic energy form. About half of this, Q — 0.25, is assumed to be radiated 
away in the external shock (forward and reverse) when deceleration starts, while 
the other half, ( ag ~ 0.25, continues as kinetic energy of the decelerating ejecta, 
and is eventually radiated away in the longer afterglow phase. Thus, ( r + (k = 1- 
Most of our results are evaluated for ( r = 1 — ( k = 0.5, and we use these nominal 
numerical fractions unless stated otherwise. These factors are used to define the 
respective luminosities, e.g., Lfc,53 = (1 — ( r )L tt5 s, etc., and we use the Q = 10 X Q X 
convention in cgs units, except, e.g. for the coasting Lorentz factor (77 = 600 rfaoo)- 

We consider a continuous outflow in which the central engine is active for 
£o,obs = 20 ti, 3 (l + z)/2 s, while the deceleration time is ~ few seconds resulting in a 
modestly thick-shell deceleration region, which for simplicity will be approximated 
as a thin shell. We assume a constant density external medium (n = 1 no cm -3 ) 
and a burst with a fiducial redshift and luminosity distance of z = 1 or = 
2 x 10 28 D28.3 cm. 

In magnetically dominated models the magnetic fields become rapidly trans- 



verse, in a striped magnetic structure (Kennel & Coroniti 1984). The Lorentz 



factor in such jets will increase more slowly than in a matter dominated case, e.g. 
after a few characteristic lengths above r the average Lorentz factor will increase 
as r m , where m depends on the geometry assumed at the lower boundary (e.g., 



McKinney & Uzdensky 2011 Narayan et al. 2010). When magnetic reconnection 



plays a role, the expected Lorentz factor behavior is 7 oc r 1 / 3 until it saturates 
( |Drenkhahn||2002| |Giannios fc Spruit||2007t |Metzger et aLpOllbt |Meszaros fc Rees 



2011), which we assume in what follows. 



Th) oc 



r i/3 if r < r sat 
const, if r sat < r < r dec , 



where r sat is the saturation radius where r — > 77 = L t /Mc 2 and r dec > r sat is 
the deceleration radius. We will concern ourselves here only with the behavior 
before and in the neighborhood of the deceleration radius, and do not consider 
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the longer term afterglow phase. A feature worth stressing is that an acceleration 
behavior such as equ. ([!]) implies a comoving volume element V oc r 2 T(r) scaling 
leading to a radius dependence of the comoving particle density n', temperature 
T", etc. which differs from that in the usual baryon-dominated 7 oc r. Also, the 
saturation radius occurs at significantly larger radii than in the matter dominated 
acceleration cases, if other relevant quantities are the same. Furthermore, the 
scattering photospheric radius (see equ. pj) occurs generally before the saturation 
radius. 



3. Prompt Photospheric Radiation 

The photospheric radius is obtained by setting the scattering optical depth 
to unity for the ejecta, i.e. t t = L r f (47impC 3 r]T)iT T r / \2Y) = 1, if the e + e _ pair 
contribution can be ignored. The effect of pairs will be that they will increase the 
photospheric radius. With the increased radius variables have to be recomputed, 
resulting in a different cutoff energy and in a different number of pairs. At the 
end of this iterative process, we find that the photospheric radius increases by a 



factor of ~ 5 for the nominal parameters (see also Bosnjak & Kumar 2012). We 



address the details of pair creation at the end of this section. 

For T = (r/Vo) 1 / 3 the photospheric radius at which Tt = 1 is 

' t 2/3 \ 3 / 5 
L r r ' a T \ 



r >* = I i^J = 6 - 5 x 10 ' 2 cm < 2 > 



for a nominal 77 = 600, where gt is the Thomson cross section. The photospheric 
radius occurs in the accelerating phase for any realistic set of parameters, for the 
magnetized dynamics described in the previous section. The saturation radius, 
where the acceleration ceases and the ejecta starts to coast with r ~ 77 ~ constant 
is, for a magnetically dominated jet, at 



sat 



r oV 3 = 2.2 x 10 15 r 0J rj 3 60O cm. (3) 



While compared to baryonic outflow models the photosphere is relatively farther 
away from the central engine, both the photosphere and the saturation radius 
remain generally smaller than the deceleration radius, 

r ^ = (i3p) = 48 * 10"4Stt - cm. (4) 
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This is the radius at which we calculate, in ^4], the interaction of the photospheric 
photons with the external shock electrons, resulting in SSC and EIC components. 
Using the scaling relations of the magnetic dynamics, the density decreases with 
radius as n' b oc V'^ 1 oc r~ 7 / 3 for r < r sat = r r] 3 and n' b oc r~ 2 for r sat < r, that 
is, n' b (r) = n^r/ro) -7 / 3 for r < r sat and n' b (r) = n' sat (r /r sat )~ 2 for r sat < r. Thus, 
the comoving baryon density at the deceleration radius is n' bdec = 5.3 x 10 4 (1 — 
Cr) 2 ^ 3 L t ^ 3 r] 60( { nj t 13 f cm 3 . 



3.1. Nonthermal Photospheric Spectrum 

At the scattering photosphere, a prompt spectrum with a high radiation 
efficiency can be produced through dissipative effects. For non-magnetic out- 



flows, such photospheres were considered by, e.g. Rees & Meszaros (2005); Pe'er 



et al. (2006); Beloborodov (2010), while for magnetic photospheres dissipation of 



magnetic energy via reconnection, turbulence dissipation or the associated semi- 
relativistic shocks was considered by, e.g. Giannios & Spruit (2007); Giannios 



(2011); McKinney & Uzdensky (2011); Meszaros & Rees (2011). This generally 



can convert a fraction ( r < 0.5 of the initial jet energy L t into prompt burst radi- 
ation, the rest remaining available as kinetic energy of the coasting ejecta beyond 
the saturation radius. The spectrum of an unmagnetized dissipative photosphere 



has been calculated by Beloborodov (2010), both analytically and numerically 



showing that a Band function spectrum is produced, with spectral indices in the 



range of those observed (see also Vurm et al. (2011) for inclusion of magnetic 
fields up to strengths eg ~ 1). A magnetic photosphere spectrum was considered 



by Thompson (1994), where reconnection leads to magnetic turbulence, resulting 
in a Band spectrum with the canonical a = 1 and /3 — 2. A magnetic dissipative 



photospheric spectrum was also calculated by Giannios & Spruit (2007), similarly 



resulting in a Band-like spectrum. Here we do not specialize to a specific dissi- 
pation mechanism, assuming that a Band function is likely to be produced at the 
photosphere. 

The value of the Lorentz factor at the photosphere is 



ph 



ph 



/ro) 1 / 3 = 87 Ll%C /5 r^%^ 5 . 



(5) 



At this photosphere, the magnetic dissipation and/or collisional effects such as 
n,p decoupling will lead to a significant non-thermal electron component. E.g. for 
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reconnection acceleration or semi-relativistic shocks associated with reconnection 



(Meszaros & Rees 2011), a typical minimum comoving electron random Lorentz 



factor 7e P /i ~ 600r r , coincidentally of similar order of magnitude as r\. The mag- 



netic field at the photosphere is B' ph = (327re,empC 2 n'fc) 1//2 rv = 1.7 x 10 6 L t ^ (1 — 
CrY^Cr 7 ^ 1C>r o 7 vloo e Bo^r G. Here n' b is the comoving baryon density evaluated 



at the photosphere, n' bph = L k / (4:7im p c 3 rp h r]T ph ) = 2.0 x 10 13 L t ^°( 1 
Cr) r o7 ^eoo cm -3 , Here as before, L k and L r are expressed as function 
( r . The observed synchrotron spectral peak of the photosphere will be 



2/5>-7/5, 



^ph,obs 
~-sy 



3ehB' 



£br 



ph 



l + Z 



3ioc 1/2 (i-c,) 1/2 r;/ 7 2 4 /2 r; 



1 + z 



keV, 
(6) 

corresponding to the Band peak. The magnetic field energy fraction at the pho- 
tosphere is still 6b ~ 1, but dissipation will lead to values €b,fs ^ 0.1 by the time 



Thompson 


( 


1994 


); 


Giannios 



& Spruit (2007) and Beloborodov (2010), we assume that the emergent spectral 



shape is of the Band form, taking for its peak value dependence on the flow pa- 
rameters the expression (J6|. As a nominal photon number low energy spectral 
index we take a ~ 1 (where N £ (e) oc e~ a ) and for the high energy spectral index 
(3 ~ 2.4. The flux in energy per energy units is then 



F E = A 



-M ife<e br 

c br J 

\ -0+1 

t) ii£>£ >»' 



(7) 



The flux is related to the luminosity via L/AitD 2 l = J F e de which yields the normal- 



ization factor A = L r /(47rD 2 L e p £ obs (l/(2-a) + l/(f3-2))) = 9.1 xlO" 9 L^Q' (1- 
C r O" 1/2 ^ L 2 8 . 3 rJ > / 7 2 e B 1 /2 r j : 3 ((l + ^)/2) erg cm" 2 s" 1 keV -1 . The eF £ peak of the non- 



thermal spectrum will arise at and at (sF £ 



e br A = 2.8x10 L t &( r D L 7& x erg cm 



e SA « 6.8 Ll%(r 1/5 (1 



At low energies a self-absorption break is ex pected. In our case, this will occur at 

Cr) 1/5 %oo 5 i y S keV (Guetta & Granot||2003a|. 



The high energy branch of the photospheric spectrum ([7]) extends up to an 
energy which is model dependent. We discuss here two scenarios. In the first 
scenario, photons above the photospheric spectral peak ^ are upscattered as a 
result of interactions with electrons associated with magnetic turbulent waves, up 
to an energy e' ~ m e c 2 in the jet frame (Thompson 1994). In our case, this 
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corresponds to an observer frame cutoff at 



e h « (A/3)m e c% h /(l + z) « 30 L\%£ V,,: 5 / ' ~ 



^600 



MeV. (8) 



This corresponds to e' ~ m e c 2 in the comoving frame, and is below the 77 cut- 
off expected from pair production against photospheric lower energy photons (see 
below), so no pairs are created. 

In the second scenario, we assume that the magnetic reconnection regions or 
the shocks associated with them are coherent over long enough times that electrons 
can be accelerated, e.g. via a Fermi mechanism, to a power law extending above 
the previous Eh to photon energies sufficiently high to be subject to 77 interactions. 
We estimate the 77 annihilation energy from requiring the pair optical depth to 
be unity against target photons of energy e at . 



-77 



(T ph m e c 2 Y/[(l + z) 2 e at ], 



(9) 



11/180^^/4^ 



(Lithwick & Sari 



2001 



Murase & Ioka 



2008), 



where r 77 ~ 

where N >£ is the number of photons with energies higher than e. This will result in 
an observer-frame cut-off energy e 77 ~ 30 — 100 MeV. If e 77 > e at is not satisfied, 
the cutoff will be at ~ r p / l m e c 2 /(l + z). Above the cut-off, the spectrum becomes 
steeper, the exact post cut-off spectral index depending on the spectrum and 
photon spatial distribution. A simple slab approximation, as discussed for GRB 



090926A (Ackermann & the Fermi collaboration 2011), results in a steepening of 



the high energy slope index by f3 — 1 above the £ 77 energy. 

These pairs also radiate in the magnetic field of the prompt emission site, and 
this will result in a low energy synchrotron component peaking at tens of electron 
volts. As we discuss in §[6j this component is a good candidate for the bright 
prompt optical emission observed in some bursts. The pair synchrotron component 



will be in the fast cooling regime, with 7±, m 
and £±, P eak ~ 3.2 eV (e.g. Toma et al. 2011a) 



3.2 eV (e.g. 

different depending on the different cases for e 77 
is 7.7 x 10 -9 erg cm" 2 s" 1 for nominal parameters. 



e 77 (l + z)/(2T ph m e c 2 ) = 3.9 
The functional dependence is 
The peak of the eF £ spectrum 



3.2. Thermal Component of the Photospheric Spectrum 



In the presence of collisional or magnetic dissipation, in addition to a non- 
thermal component one expects also a thermal component, whose luminosity 
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Lth = CthLt should have a quasi-blackbody spectrum. In our case, this peaks in 
the soft X-rays (Meszaros & Rees 2011), due to the different temperature scaling 
with radius for the magnetic dynamics, comprising a fraction Cth of the luminosity 
of the photosphere. This thermal component can be calculated from the initial 



T 



1/4 



2.1 L\ 



1/4 



,53 r 0,7 



l/2 r -l/2 
1 r 



MeV at the initial radius tq = 10 7 ro.7 cm. 



Above r the magnetically dominated jet dynamics T oc r 1//3 implies a comov- 
ing volume is V oc r 2 T oc r 7 / 3 in the acceleration regime. Thus, the temperature 
will decrease more gradually with radius than in the baryon-dominated case (where 



T' oct 1 ), as T' oc p 



n 



1 oc p' 1 / 3 oc V- 



1/3 



oc r 



7//9 , or T'{r) = T (r/r ) 7 ^ 9 , where 



7 = 4/3 is the adiabatic exponent for a relativistic gas. Thus, at the photosphere 
r(r P = 27L-rCr 1/15 ,r%P°r r - 1 ' 2 (¥)" 1 keV. 



The corresponding thermal luminosity is L t h = 6.5 x 10 49 erg/s, which is low 
for producing the prompt emission, as well as being too soft. This blackbody 
component from the photosphere thus peaks in the soft X-rays at (sF Et BB) peak ~ 



9.5 x 10- 9 L t U / 3 15 C; 



■4/15 4/15 4/15 
'/600 '0,7 



erg cm s . This thermal component is similar 



to the one found by Page et al. (2011). 



4. External Shock Radiation Spectral Components 

4.1. Forward Shock (FS) Synchrotron 

The forward shock develops at the deceleration radius, where the jet has 
plowed up an amount of external mass roughly equal to 1/r] times of the ejecta 
mass. The electrons in the shock will be accelerated into a relativistic energy 
distribution, and will undergo cooling through synchrotron emission and by in- 
verse Compton scattering off external (photospheric) and their own (synchrotron) 
photons. At this deceleration radius the Lorentz factor has roughly halved from 
its coasting value. The time of the deceleration is td ec = r dec/(2^ 2 c)(l + z) ~ 
4-4 L)%{1 - CrY^n-^rj-^il + z)/2 s. 

While at the photosphere the magnetic field parameter e p ^ is close to unity 
for a magnetically dominated model, after magnetic dissipation ceases one expects 
this parameter to be much less, the outflow becoming essentially baryon dominated 



(e.g., 
0.1. 



Zhang & Yan 2011). Here we assume that at the deceleration radius e^ s ^ 
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The magnetic field in the forward shock is then B' FS = (327im p c 2 eBT 2 7S n) 1 ^ 2 « 
74 ^600 e s 2 -i n o^ 2 The cooling Lorentz factor is 7 C = Qirm e c/ ((l+Y)aTB FS Ttde C ) = 
3/(8(1 + yyrK/mJeB nr^T?) = 98 (1 - Cr)~ 1/3 L;^ 3 h^e^n, 2/ %^{l + 
Y)" 1 , where Y = (—1 + a/1 + Ae e /eB,Fs)/2 ~ 0.1 is the Compton parameter. The 
minimal (injection) Lorentz factor is 

7m = ej^-- — \y = 3100 e e .277600^,2.4, (10) 
m e p — 1 

where g p> 2A is the fraction (p— 2) / (p— 1) normalized to p = 2.4. Thus, the electrons 
responsible for the synchrotron radiation are in the fast cooling regime, and their 
distribution is given by dN e /d^ oc 7~ 2 if 7 C < 7 < 7 m and dN e /d^ oc 7~ p ~ 1 if 

7 > 7m- 

The observer frame energy of photons radiated by electrons with a random 
Lorentz factor 7 m in the forward shock is 

_ 3heB' FS 2 _ 1/2 2 1/2 4 2 f l + z \ 1 , s 

£ m — c 7m _|_ ^ — 6 -° e B ,-l e e ,-2 n 7600^,2.4 1 2 J ' ^ ' 

while for electrons with a random Lorentz factor 7 C the observed photon energy is 

e c = 3.7 L- 5 f(l - W-^e^no^Cl + >T 2 (^) eV. (12) 

The peak flux density of the FS synchrotron spectrum occurs, for fast cooling, at 
e c , and is given by (e.g. Sari et al. 1998): 

f fs (r x _ 4irr 3 dec nm e c 2 a T r)B FS _ , , 1/21/2 2 Tv n ^ 

^max{ £ c) - g 127Tg D 2 ~~ L t,S3{i- - (,r)tl.3^B,~l n U L,28.Z J Y' \ 16 ) 

The maximum Lorentz factor attainable by the electrons is calculated by equating 
the acceleration timescale in the shock to the radiation timescale, which gives = 
(3q e / gM&TB' FS ) l l 2 ~ 5.4 x 10 6 Vmi 2 ^-! 71 ^ 1 ■> where gu is a numerical factor of 
order unity. The corresponding photon energy iseps,M — 3hq e B' FS / (4nm e c)-y 2 I ri/(l + 
z) = 11.2 ?76oo((l + -2)/2) _1 GeV. We note that the maximum synchrotron energy 
only depends on the terminal Lorentz factor. 

The scattering optical depth in the FS is given by 

N e a T 4nr dec na T 8 1/3 0/3.1/3 -2/3 2/3 n ,x 

where N P is the number of electrons in the forward shock. 
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4.2. Reverse Shock (RS) Synchrotron 

The reverse shock, if it develops, becomes strongest at the deceleration radius. 
A contact discontinuity (CD) separates the FS from the RS, and in the frame of 
the CD the RS travels backwards. The pressure is the same in both sides of the 
CD, and consequently the magnetic energy density will be the same as well, under 
the usual assumption that the magnetic field is turbulently generated and is 
taken to be the same in both sides. 



The characteristic frequencies for the RS will be: e 



RS 



-FS 



ID" 5 eZA-A^oodhA^T 1 keV and s** = (e.g. 



1/2 9 



HI 



1.1 X 



Sari & Meszaros 



2000). The RS cooling frequency can differ by about 10%, but we will ignore this 
difference. Generally the RS peak flux will be a factor V higher than in the FS. In 
our case the radiative regime changes from the FS to RS: while the FS electrons 
are in fast cooling, the reverse shock is in the slow cooling case. We calculate the 
peak flux from 



f rs t 



N P P n 



Lkto m e c 2 aTf]B 1 



4tiD 2 l 



3<? e 



92 Li )5 3(l-Cr)tl.3f/600eB 



1/2 



l u -^£,28.3 J J- 



(15) 

It is unclear, in an initially magnetically dominated outflow, whether the RS will 



develop or not, depending on various assumptions ( 


Zhang & Kobayashi 


2005 


Giannios et al. 


2008; 


Zhang & Yan 


2011 


Narayan et al. 


2011 


). In the simplest 



cases the Alfvenic 7^ ~ a/1 + a Lorentz factor will be high and magnetic waves 
can carry information from the CD to the start of the ejecta, suppressing the RS. 
Here a ~ F £>r /F £ ^ is the magnetization parameter giving the ratio of the Poynting 
(magnetic) flux to the kinetic flux. However, the result depends on the amount 
of dissipated magnetic energy at radii before reaching the deceleration radius. If 
there is significant portion of magnetic energy compared to the kinetic energy there 
may be no reverse shock, while in the contrary case reverse shocks may form. 

The optical depth of the RS region is important for calculating the flux of the 
inverse Compton radiation. We calculate the optical depth of the RS at r^c from 



trs 



N e a T 



4tiR 2 r)m p c 2 <i'Kr 2 :ec 



6.4xlO^(l-C) 1 /M / ^o / o 3 -o /3 (16) 



where N e is the number of electrons in the reverse shock. Note that at r^ ec the 
optical depth of the RS is 77 times the optical depth of the FS as expected. 
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4.3. Forward Shock Self-Compton (FS-SSC) 

At the deceleration radius the forward shock (and the reverse shock, if present) 
accelerated electrons will cool also by inverse Compton interactions with their 
own synchrotron photons. The forward shock at the deceleration is in the syn- 
chrotron fast cooling regime, for the parameters considered here, and the Comp- 
ton parameter is Y = (-1 + y/l + A-qssc^-el Vb) /2 « e e /e B = 0.1, where rjssc = 
min((7 c /7 m ) 2_p , 1) and rj SS c = 1 is valid for the fast cooling case. 

The SSC component of the FS will also be in the fast cooling radiative regime, 
and the minimal and the cooling energies will be given by £^f ~ ^1rn £ m = 
75 ^ 2 -i<- 2 nf% 6 ooSi4((l + z)/l>T l GeV and ef c « 2 7c 2 e c = 70 L~ 5 4 3 /3 (l - 
Cr)- A/ H^%]iW 3/6 vS^ + >T 4 ((1 + ^/S)- 1 keV respectively. 

The amplitude of the FS SSC component is determined from: (e ssc F E )^^ = 

Ys^FfJ = 3.7 x 10-% 2/ 5 3 3 (l - C) 2/3 ^ / 3 3 % 8 o / o 3 n 1 /3e^ 2 < 2 . 4 ^ 2 2 8.3 erg cm' 2 s~\ 
The higher energy electrons can reach the Klein-Nishina (KN) regime. This re- 
sults in a break in the spectrum at Sjf^ c = 7 c m e c 2 r]/(l + z) = 15 L^( 3 (l — 

Cr)" 1/3 *r.3 /3e B-i w o 1/3 %oo( 1 + Yssc)' 1 GeV. Details of the KN break energy are 
in Appendix [B} 



4.4. EIC scattering of photospheric nonthermal photons on external 

forward shock electrons (FS-EIC) 

The electrons in the forward shock will also lose energy by external inverse 
Compton (EIC) as they upscatter the photons from the prompt emission (photo- 
spheric) region. The minimum and cooling frequencies of this component will be 
e E J C « 2 7 ^ br = 6.2 Cr" 1/2 (l-C,) 1/2 4(o<-2% 2 oo^2.4^ /2 r 3 ((l + ^)/2)- 1 TeV and 
e*'° ~ 2^ c e br = 6.0L- 5 f C r 1/2 (l-C,)-^-f ^-i^W^of (^-^^((^ 
z)/2) 1 GeV, where we differentiate between prompt and FS magnetic parameter 
by using e B an d cb-i respectively. 

In the absence of KN suppression, the flux of the FS EIC emission would be 



of the form (Murase et al. 2011): 



{e 1 a ife < £ E ic,c 
e- 1 / 2 ife EIC , c <e<e EIC , m (17) 
e~ p/2 iie E ic,ra<e 
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The KN effects however would introduce breaks in the high-energy part of the 
spectrum. Details of the derivation of the KN frequencies are in Appendix [Bj For 
a large part of the parameter space, £Eic,m will be in the KN regime, and the 
eF £ peak will be at the KN break frequency e^ 1 ^ ~ ^7 c m e c 2 /(l + z) = 15 (1 — 
Cr)- 1/3 ^ /3 ^J /3 eiJ-i^ 2/3 ^(l + Y)~\(l + z)/2Y l GeV. A rough estimate of 
the peak of the EIC emission at the peak is 

{eF £ y«% ~ e br N £iP r FS e^ = (18) 
7.2 x 10- 10 L t , 53 C, 3/2 (l - Cr)- 1/2 6^ /2 ^-ir r - 3 (l + Y)^Dll %2 erg cm" 2 s^W) 

where N £jP is the photon number spectrum of the prompt emission, evaluated here 
for a — 1, /3 = 2.4 and for e ~ e^. To account for the anisotropy of the emitted 



radiation in the forward shock frame, we multiply our flux by a factor of 0.5 (Fan 



& Piran 2006). 



4.5. EIC scattering of photospheric non-thermal photons on external 

reverse shock electrons (RS-EIC) 

If the reverse shock develops, there will be an external inverse Compton 
component from the reverse shock electrons scattering prompt photons as well. 
This differs from the forward shock EIC because of the larger optical depth and 
the lower energy of the electrons in the RS, and because the RS will be in the 
slow cooling phase. The characteristic frequencies are e^f EIC ~ ^lRSm £ br — 
ire e ,_ 2 ^ 2 .4ro, 7 e^ /2 o r3((l+z)/2)- 1 MeVandef^ « 2 7 ^ c e br = 6.0 L~%\r 1/2 (1- 

Q- l/ %T 4>^W% 6 ~oo /3 (l + Y)- 2 r^ 7 /2 T%(l + zJ/2)- 1 GeV. The RS elec- 
trons, being in the slow cooling phase, will upscatter the prompt emission into a 
spectrum of the following shape: 

e l ~ a if e < e* SEIC 

pRSEIC ^ ^ £ (l- P )/2 i{£ RSEIC < £ < £ RSEIC (2Q) 
£ -p/2 ii£ RSEIC < £ 

The peak of the emission will be at the Klein Nishina cutoff frequency for 
this component, which is e^f EIC ~ rj^^fm e c 2 / '(1 + 2) ~ 0.8 ?760oe e ,-25'p,2.4((l + 
z)/2) _1 GeV. This introduces a spectral break and the photon index above this 
energy will be ~ — (a — p — 2). 
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The peak flux of the emission, considering a weakening by 0.5 due to anisotropy, 

is 

eF$BBic - e br N, p r RS ef N ^ = 2.3xl(T 8 ^^(l-^)"^^ 

(21) 

Note that e KN ' m falls in the low-energy part of the Fermi LAT range. Thus, if 
the RS-EIC were dominant, the photon index expected would be — (a— p— 2) rj 3.4. 
In bursts with extra high-energy components, such soft photon indices (or even 
softer) are indeed observed. For other model parameters, we found that while 
the RS-EIC is the dominating component at its peak, other inverse Compton 
radiation components such as FS-EIC of FS-SSC can make the spectrum harder 
by contributing at a few x 10 GeV s. 



5. Numerical Results and Model Parameter Variations 

We calculated a number of model spectra based on the considerations of the 
previous sections. The initial set of nominal parameters used is L t = 10 53 erg/s, 
t = 20 s, Cr = 0.5, Cfc = 0.5, rc=l cm -3 , 7/ = 600, e B , P r = 1, e B ,FS = £b,rs = 0.1 , 
e e ,F5 = ee,RS = 0.01, r = 10 7 cm, z = 1 , D L « 2 x 10 28 cm. 

Here ( r = 0.5 is the nominal fraction of the initial L t assumed to be radiated 
at the photosphere in non-thermal and thermal radiation, Ck = 0.5 is the energy 
radiated by the external shock at the deceleration radius Also, for simplicity, when 
a reverse shock forms, we assume the same values for e# and e e in both the forward 
and in the reverse shocks. We have explored also the effects of departures from 
these various nominal parameters. 

Models with Negligible Pair Formation - A first set of models was calculated as- 
suming that the first scenario of £3.1, where the photospheric spectrum of eq. (j7l) 



cuts off around m e c 2 in the jet frame, or around ~ 50 MeV in the observer frame, 
with pair formation being negligible in both the photosphere and in the external 
shock. The resulting generic observer-frame spectrum consists of a Band func- 
tion spectral component peaking at sub-MeV energies and extending up to ~ 50 
MeV, with a second component at GeV energies due to a combination of EIC and 
SSC by external forward and reverse shock electrons of photospheric photons and 
their own synchrotron photons. This second GeV component has typically a total 
fluence which is ^ 0.1 that of the MeV Band function, and for some parameter 
ranges it stands out from the first component, while for others it merges more 
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or less smoothly with the high energy branch of the first component. Whether it 
stands out or not depends on the external density, the external shock parameters, 
and on whether a reverse shock forms or not. 

10 6 i 1 




e [MeV] 



Fig. 1. — A model without pair formation, L t = 5 x 10 52 erg/s, t = 20 s, ( r = 
0.6, G = 0.4, n = 100 cm" 3 , r) = 400,e BjPr = 0.9,e B , F5 = 1 x 10~ 2 , e e , FS = 
2 x 10~ 2 , ro = 10 7 cm, z = 1, (3 = 2.4, p = 2.4. The black dashed line is the prompt 
synchrotron emission, black thin continuous line is the prompt thermal component 
(marked BB), the thick black line is the RS-EIC, the gray, thick, dotted line is the 
forward shock synchrotron part (FS), the gray, dashed line is the forward shock 
external inverse Compton, the black dash-dotted is the FS-SSC component, the 
gray, dashed line is the RS-SSC and the gray, dash-dotted is the RS synchrotron 
component. The thick gray continuous line is the sum of the components (the 
upper one is with and the lower one without the RS contributions). 

In Figure[T]we show one of the cases where the second, GeV component, stands 
out from the first Band component. The parameters of the external forward shock 
are e# = 0.01 and e e = 0.02, and the dissipative photosphere produces a Band 
function peaking at < 1 MeV in the GBM range. 
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In the LAT range, if a reverse shock is present, we can see a clear "bump" 
in the spectrum around 1 GeV (the upper thick curve). The prompt emission 
and the bump up to ~ 3 GeV in Figure [T] is strikingly similar to the spectrum of 
GRB 090926A ( |Ackermann &: the Fermi collaboration||2011[ their fi ffure 5.), which 
showed an extra power law with a cutoff besides the Band component. Here the 
extra power law would be the rising part of the RS-EIC (the model photon index 
{—(—p — l)/2 ~ 1.7 coincides precisely with the measured one for p = 2.4) and 
the cutoff is the part after the RSEIC peak. The ratio between the GBM and 



LAT ffuences is ~ 10 1 , an average ratio in the observed bursts (Omodei & the 



Fermi LAT collaboration] 2011 Pelassa 2011). At optical wavelengths the reverse 



shock produces a high flux (~ 3 x 10 8 erg cm 2 s 1 , or ~ 7) which could be 
responsible for the very bright optical flashes observed in some bursts. 

However, if the reverse shock is weak or missing, in the same figure [T]the upper 
thick bump at ~ 1 GeV is absent, and is replaced by the lower thick line. One has 
the same GBM prompt emission and a hard component emerging at < 10 GeV. 
This component would also appear as an extra power-law component, but in this 
case, it is weaker by one order of magnitude than the prompt emission. For the low 
photon number detected in this energy range, such a component might have low 
significance in a fit. Also the optical flux is much lower, 3 x 10" 10 erg cm -2 s" 1 , 
or rrifi ~ 11.8. We note that the combination of parameters used for obtaining the 
single component and two component spectra in this figure are not unique. 

In Figure [2] we show one of the cases where the overall spectrum appears as 
a single Band component, due to the second GeV component approximating an 
extension of the first, MeV Band component. The parameters of the external shock 
are in this case e# = 0.02, e e = 0.005 for both the forward and reverse shock, with 
r) = 400 and n = 30 cm" 3 , other parameters being nominal. It is seen that the 
combination of the external reverse shock EIC and the forward shock SSC plus 
EIC provides a roughly continuous slope connecting the MeV Band spectrum to 
the high energy branch (the upper thick curve from ~ 0.1 GeV onwards). The 
reverse shock EIC hump fills in the through left by the forward shock SSC. The 
wiggles from the two humps at 100 MeV and 10 GeV would be smoothed out in 
typical spectral fits with low photon number statistics in this range encountered 
in Fermi LAT bursts. The matching of the flux levels between the photospheric 
spectrum high energy end and the external RS peak is controlled largely by the 
external density n, whose value here is in the usual range. 



s [MeV] 



Fig. 2. — A different model, also without pair formation, and parameters L t = 
10 53 erg/s,t = 20 s, ( r = 0.5, (k = 0.5, n = 30 cm" 3 ,?? = 400,e B , pr = l,e B ,FS = 
e B ,RS = 2 x 10- 2 ,e e , FS = e e>RS = 5 x lO" 3 ,^ = 10 7 cm, z = 1,(3 = 2.5, p = 2.4. 
The black dashed line is the prompt synchrotron emission, black thin continuous 
line is the prompt thermal component (marked BB), the thick black line is the RS- 
EIC, the gray, thick, dotted line is the forward shock synchrotron part (FS), the 
gray, dashed line is the forward shock external inverse Compton, the black dash- 
dotted is the FS-SSC component and the gray, dash-dotted is the RS synchrotron 
component. The thick gray continuous line is the sum of the components (the 
upper one is with and the lower one without the RS contributions). 
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Assuming that a reverse shock is absent or negligibly weak (i.e. ignoring the 
upper thick gray continuous line in fig. [2] peaking at > 100 MeV), one sees that 
now a second hard component appears above 100 MeV and peaks at 3 GeV, being 
in this case mainly due to the forward shock SSC, with a contribution from the 
EIC of photospheric photons. Other components, while present, are here sub- 
dominant. The secondary peak (also in the absence of the RS-EIC component in 
fig. |2]), though reminiscent of the extra power-law component observed in GRB 
090902B or 090926A, is ~ 2 orders of magnitude fainter than the main peak, while 
in the observed cases this ratio, when reported, is at most ~ 1.5. The optical flux 
predictions are w,r ~ 8.3 and rriR ~ 10.5 in the presence or absence of the reverse 
shock. Again, this is not the sole combination of parameters which produces an 
approximate single Band function extending to GeV energies. 

In the above models where pair formation is not expected, the qualitative 
effects of increasing the terminal Lorentz factor 77 consists in a strengthening of 
the RS-EIC component. Even though at low rj we expect a lower cutoff for the 
prompt emission and a more prominent inverse Compton component, this is not 
the case. Below 77 ~ 300 the main high-energy component is the sum of the FS- 
SSC, the FS-EIC and the prompt SSC at ~ 10 GeV which form a distinct peak. 
Increasing the density n makes the RS-EIC component to be more detectable 
from n ~ 10 cm 3 up to n ~ 500 cm 3 . Outside these parameters the FS-SSC 
and the FS-EIC are the same magnitude or dominating the RS-EIC component 
and a break in the spectrum is more visible. The shape of the high energy part 
of the spectrum has only a weak dependence on the value of e#. By decreasing 
the prompt magnetic parameter, the peak energy becomes lower and the external 
components could become more prominent. 

Models with Pair Formation- A second set of models was calculated assuming 
that the Band spectrum from the photosphere extends to sufficiently high energies 
that pair formation from 77 interactions is expected (see In these cases the 
primary Band component extends up to an energy given by eq. ([£]), and there 
is a secondary synchrotron component from the 77 pairs, both the primary and 
secondary photons leading to separate SSC and EIC components from upscattering 
in the photosphere and in the forward and reverse external shock, some components 
being more important than others. Figure [3] presents one of the cases involving 
pair production. The overall behavior is not too different from that of the models 
with no pair production: the RS-EIC component is approximately smoothly joined 
to the prompt at ~ 100 MeV. The RS-EIC in turn joins smoothly to the other 
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Fig. 3. — Model with pair formation with L t = 10 53 erg/s, t = 20 s, ( r = 0.6, (k = 
0.4, n = 10 crrr 3 ,?? = 600, e B , P r = i,e B ,FS = £b,rs = 1 x 10~ 2 ,e etFS = e etRS = 
1 x 10~ 2 , r = 10 7 cm, z = 1,(3 = 2.4, p = 2.4. The black dashed line is the prompt 
synchrotron emission, black thin continuous line is the prompt thermal component 
(marked BB), the thick black line is the RS-EIC, the gray, thick, dotted line is the 
forward shock synchrotron part (FS), the gray, dashed line is the forward shock 
external inverse Compton, the black dash-dotted is the FS-SSC component, the 
gray, dash-dotted is the RS synchrotron component and the thick dash-dotted 
is the pair synchrotron contribution. The thick gray continuous line is the sum 
of the components (the upper one is with and the lower one without the RS 
contributions) . 
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higher energy components at ~ 10 GeV. If a RS does not develop, the only missing 
component at high energies will be the RS-EIC (thick black line) and again we 
get a bump in the spectra at ~ 10 GeV energies. The optical flux is only about 
1.3 magnitudes fainter in the absence than in the presence of the reverse shock, 
771r ~ 9.3 and ~ 8 for the two cases. 

In Figure [4] we show another case where pair formation occurs, for a different 
choice of parameters. In this case the second component bump would be present 
irrespective of whether the reverse shock is present or absent, since both shocks 
result in a similar bump. The FS synchrotron and the FS-EIC contribute signif- 
icant flux to the bump. In this case, the optical flux with the reverse shock is 
larger, ~ 5.5, and much fainter in the absence of the reverse shock, m# ~ 10. 

Low Energy Power Law Extensions - A notable feature of Figs. [T] through [4] is 
that the forward shock synchrotron radiation extends into the optical range. In 



Asano et al. (2010), such a power law extending into the optical was obtained for 
a hadronic cascade GRB model. Here it arises in a purely leptonic mode. For 
reasonable parameters this component falls between the cooling and the charac- 
teristic synchrotron frequencies. One sees that, for bright Fermi LAT bursts, in 
the absence of pair formation the forward shock synchrotron can produce a prompt 
optical flash of > 11 — 12, even in the absence of a reverse shock. This vi- 
sual flux scales with the prompt photon luminosity L r . A conversion of the flux 
units to R-magnitudes is 10 -8 erg cm -2 s _1 — > ~ 8.0, and 10 -10 erg cm -2 
s _1 — > itlr ~ 13.0, with the usual five magnitudes interval per each factor 100 
increase in flux. In the cases where pair formation occurs, the pairs contribute 
an additional synchrotron component, which in the optical range predicts an even 
brighter prompt flash. This is seen in Fig. |3j with the source at z — 1. In both 
types of models (with or without pair formation), in the cases where a reverse shock 
is present the optical flashes can be significantly brighter, in the range ttir ~ 6 — 8 
at the same redshift. 



6. Discussion 

We have addressed the high-energy spectral properties of the bursts observed 
with Fermi LAT , using a magnetically dominated outflow model where the prompt 
MeV emission arises in the photosphere, and high energy components arise from 
inverse Compton scattering by both photospheric and external shock electrons. 
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Fig. 4. — Another model with pair formation, and parameters L t = 10 53 erg/s,t = 
20 s, Cr = 0.2, C fc = 0.8, n = l cm~ 3 , rj = 600, e BjPr = 1, € B ,fs = £b,rs = 0.1, e ejjF s = 
e e ,R5 = 1 x 10" 2 ,r = 10 7 cm, z = = 2.4, p = 2.4. The black dashed line is 
the prompt synchrotron emission, black thin continuous line is the prompt thermal 
component (marked BB), the thick black line is the RS-EIC, the gray, thick, dotted 
line is the forward shock synchrotron part (FS), the gray, dashed line is the forward 
shock external inverse Compton, the black dash-dotted is the FS-SSC component, 
the gray, dash-dotted is the RS synchrotron component and the thick dash-dotted 
is the pair synchrotron contribution. The thick gray continuous line is the sum 
of the components (the upper one is with and the lower one without the RS 
contributions) . 
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We have investigated circumstances under which a single Band function appears 
to extend to the highest energies detected by the Fermi LAT , and where a second 
high energy component shows itself above the MeV range Band spectrum. We have 
also addressed, in the same context, the production of prompt optical flashes. 

We find that, qualitatively, these types of models are able to explain the 
observed diversity of GBM and LAT prompt emission spectra, without invoking 
an internal shock emission at intermediate radii. The presence or relative strength 
of a reverse shock plays a role in determining the spectral shape and the basic 



dichotomy in the latter. As discussed by, e.g. Narayan et al. (2011 ); Mimica & Aloy 



(2011); Giannios et al. (2008), for modest magnetization a of the outflow a reverse 



shock may be expected; and while for initially magnetically dominated outflows a 
(or es) is still large at the photosphere, at larger radii it can be expected to drop 
sufficiently low to make reverse shocks possible, although uncertainties remain 
concerning the threshold magnetization and the reverse shock strength. Thus, 
we assumed that within the normal range of parameters, reverse shocks may be 
important in some bursts and not in others. When a reverse shock is absent, the 
forward shock results in a distinct high energy (GeV) spectral component, whose 
fluence is 1-2 orders of magnitude below that of the Band component. In these 
cases, the effective high energy slope /3 of the Band component appears rather 



soft, in agreement with data discussed in, e.g. Zhang et al. (2011). On the other 



hand, when a reverse shock is present, its inverse Compton radiation can result in 
a larger fluence extra high energy component, 0.5-1 orders of magnitude below the 
Band fluence; or it can result in a smooth continuation of the Band component, 
which can mimic a single Band high energy component of relatively hard effective 
slope (3. The relatively small number of photons observed in this band could result 
in fits where the slight wiggles in the theoretical spectra are largely ironed out. 

The fact that the LAT emission is contributed by the external shock im- 
plies that it will appear with an intrinsic delay of order t^i ~ r<fe c (l + z)/cr] 2 w 
4.4 ^53(1 — CrO^^o^o %oo (1 + s ) °f order a few seconds relative to the 
photospheric MeV component. Such a delay was indicated also in some early work 
on LAT spectra (e.g., Ghisellini et al.||2009 Kumar fc Barniol Duran||2009 ), where, 
however, the LAT emission was attributed to the forward shock synchrotron radi- 
ation. In our case, it is the inverse Compton components of the forward or reverse 
shock which dominate the LAT emission. 



The implications of our model for constraints on the bulk Lorentz factor are 
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much less stringent than in one- zone models where the GBM and LAT emission 
are assumed to arise in the same region. E.g., in|Abdo fc the Fermi collaboration 



(2009); Ackermann & the Fermi collaboration (2011); Abdo et al. (2009), such 



analyses indicated Lorentz factors 7 ~ 800 — 1000 or higher. However (Zou et al. 



2011 Pe'er 2011) in generic two-zone models the Lorentz factor need not be so 



large. Specifically, in our model, which is a two-zone model in which the high 
energy photons arise in the external shock, the compactness parameter in the 
latter is low, and the spectra can be reproduced with terminal Lorentz factors 

r ~ n ~ 300 - 600. 



As discussed in £3.1, the magnetized dissipative photosphere can produce a 



Band-like non-thermal spectrum resembling the observations. In our magnetized 
models, in the absence of pair formation this component cuts off above ~ 50 MeV, 
or in the presence of pair formation, it steepens by one power law index above 



100 MeV. Recently Kocevski & The Fermi Collaboration (2012) analysed the 



paucity of GRBs measured by LAT. They claim that nearly half of the bursts 
detected by GBM which were in the LAT field of view required a break under 
0.1 GeV to explain the nondetection by LAT. In the framework of this model 
these results mean that there is indeed a cutoff at ~ 50 — 100 MeV. This is either 
due to a softening because of pair creation or a cutoff according to the magnetic 
acceleration mechanism. 



Recently Yonetoku et al. (2011) reported a polarization measurement in the 



prompt emission from GRB 100826A. This could be an indication of a magnetically 



dominated photosphere (Waxman 2003 Nakar et al. 2003), although polarization 



might also be expected from processes not requiring strong magnetic fields (e.g., 
Lazzati et al. 2004). For a magnetic jet, the transverse field components will 



dominate in the emission region, and while for an observer line of sight along the 
jet axis the polarization could average itself out, for the larger probability off-axis 
viewing directions a net polarization could be expected. 

The magnetic photosphere models also predict a weaker thermal component 

through 111 



peaking at a few keV (£3.2), plotted in Figs. [T] through |lj Such a thermal 
component has been reported in Guiriec et al.| ( |2011 ) for GRB 100724B. In our 
models such a component appears at approximately the right energies, its fluence 
generally being low compared to the nonthermal components. However, only for 
a relatively small range of parameters would it appear possible to detect it. One 
problem is that it can be conflated with the contribution of the forward shock 
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synchrotron, e.g. as the bump around < 10 keV in Fig. [3j 

Another interesting component is the optical band extension of the external 
shock synchrotron spectra (see Figs, [I] through [1]). It is seen that just the external 
forward shock by itself already can produce optical flashes of ttir > 12 (Figs. [l}[2]), 
while if pair formation occurs in the photosphere, the cooled pairs there can lead 
to flashes of m R ~ 9.3 (Figs. [3j [4]) or even brighter for suitable sets of parameters. 
On the other hand, when a reverse shock is present, its synchrotron component 
naturally produces a bright prompt optical flash, as known for quite a while (e.g., 



Meszaros & Rees 1993, 1997). Here, in addition, we have considered also the IC 
components of the reverse shock, and the effect of a photospheric EIC component as 
well. For reasonable parameters, the flux can be close to few xlO erg cm -2 s _1 in 
the optical band, which translates to wir brighter than 7 at z = 1. Generically, the 
relative scarcity of observed optical flashes may be attributed to the fact that their 
brightness scales roughly the same way as the prompt GeV luminosity ((ei^^ak (X 
L? /3 oc L° t 8 } (eF E )l r °™ pt oc L t and (eF e )^ EIC oc L 4 t /3 ), as well as to the fact that 
reverse shocks may be rare in magnetically dominated outflows. 

For the naked eye GRB 080319, the roughly similar behavior of the optical 



and 7-ray light curves can be used to argue for a common origin of both (Racusin 



et al. 2008). However, the optical light curve is not sampled as well as the 7-ray 
light curves, and at least in the prompt phase shows temporal structures (peaks) 
comparable in duration to the deceleration time which could be compatible with 



a reverse shock origin, although the RS origin was disfavoured by Racusin et al. 



(2008) . On the other hand, a fast variability of the prompt optical flash might be 
suggestive of an origin in the same region as the prompt MeV emission, which might 
be attributed to the photospheric cooled pair synchrotron component. However, 
only very rare parameter combinations could push the optical flux of the pair 



synchrotron component up to m# ~ 5 as in the naked eye burst (Racusin et al. 



2008). E.g. in our model a combination of r\ ~ 1000, (3 w 2.1 and ( r ~ 0.99 would 



approach such brightnesses. 

We note that this model may be applicable both to long and short bursts, 
since a magnetized photosphere and an external shock would be expected in both 
cases. The relatively shorter GeV-MeV delays in the short burst cases could be 
understood in terms of a closer-in deceleration or a larger Lorentz factor, the latter 
being suggested also by their harder MeV spectra. 

Finally, we point out that the predicted inverse Compton components in these 
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models extend into the TeV range, for a range of parameters. The photon with the 



highest energy detected by LAT from a GRB had an energy of ~ 33 GeV (Abdo 



fc the Fermi collaboration| [2009). However, the spectral features above this energy 



are in the range of ground-based Cherenkov telescopes, including also HAWC and 
the future CTA, providing potentially interesting targets for such detectors. 



We acknowledge NASA NNX09AL40G, NSF PHY-0757155 and OTKA grant 
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A. Other Radiation Components 

Other electron inverse Compton components, besides those already discussed, 
may be present if a reverse shock develops. In this case one would expect scattering 
of reverse shock photons on the forward shock electrons, and FS photons on the 



RS electrons. Both of these components have a flux density (He et al. 2011) 
tfsF^L ~ ~ 10_6 J y peaking at e » 0.1 keV with a flux of {eF £ ) peixk « 

2 x 10~ 13 erg cm -2 s -1 , which makes them negligible compared to other (e.g. FS 
or RS ) components. 

The external inverse Compton radiation of the pair synchrotron photons on 
the FS and the RS will also give a flux of the order 5 x 10~ 10 erg cm -2 s _1 at 
80 MeV and 0.4 MeV respectively. 

The prompt thermal (blackbody) component will also be upscattered at the 



FS as well at the RS (Ando & Meszaros 2008). Both components consist of a set 



of smoothly joined power-laws and a cutoff at high energies. The BB-FS-EIC has 
a peak of 9 x 10 -12 erg cm' 2 s _1 at ~ 30 GeV while the BB-RS-EIC has a peak 
of 6 x 10 -12 erg cm -2 s -1 at ~ 4 GeV. 

The synchrotron self Compton component of the reverse shock peaks (in eF £ ) 
& t £ cRS ~ < ^l'cRs e c,RS ~ 35 keV. The Compton Y parameter, which gives the 
luminosity ratio of the SSC to the synchrotron component is = (—1 + 



y/1 +4r ] f t s s c e e /e B )/2 « 0.031 (Sari fc Esin 2001 ), where = min((7 c RS /7^f ) 2 ~p, 1) 



In this slow cooling regime, ?7^f c < 1 is valid. The amplitude of the RS-SSC is 
obtained from (eF £ ) p s e ^ = Y^ c ef s F e R c s « 3.5 x 10" 10 erg cm" 2 s _1 . While this 
component can have fluxes of the order of 10~ 8 erg cm' 2 s _1 , it is dominated by 
other components. 

The prompt emission will generate an SSC component. The peak energy of 
this component 2>Je,ph £ br ~ 0.2 TeV falls in the deep KN regime and its eF s peak 
will be less than 10 -9 erg cm -2 s _1 . 



B. Klein-Nishina break for SSC and EIC components 

Klein-Nishina effects are potentially important when assessing high-energy, 
inverse Compton components. In our study we investigate the KN break for SSC 
and EIC components of both the FS and RS with the prompt emission. We 
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calculate the KN break energy based on Guetta & Granot (2003b). The KN break 
occurs at the solution of the e = •y max (e)m e c 2 , where "f ma x is a function of e (for 
the detailed expression of 7 maa: (£) see Guetta & Granot 2003b). Depending on 
the radiative regime and the position of e^ N with respect to the characteristic IC 
frequencies, we have four cases for each regime (IC here stands for either EIC and 
SSC emission). In the fast cooling case (valid here for FS-EIC and FS-SSC) we 
need to solve: 



tic 

KN 



{m e c 2 ) 2 /e' c 

(m e c 2 ) 2 /e' m 
m e c 2 lM 



if 
if 
if 
if 



W < £ ,IC 

7 c fc c ^ fc KN 



< 7 2 e' 



me 

ry? pi p' IC <T >y 2 p' ( 

lm c c - c KN - lm c m\ 
1m £ 'm < £ KN < 1 M £ 'm 
lM t m - b KN 



JIC\ 



(Bl) 



for e' KN . Above the break the spectrum will change to F e oc or in some 

cases to F £ oc e~( 2 ~ a ) in the fast cooling, or to F e oc in slow cooling regime. 

In the slow cooling regime (used here for RS-SSC and RS-EIC) we can obtain 
the break energy by swapping the roles of the cooling (c) and characteristic (m) 
quantities in the above equation. 
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